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Abstract. We propose the realization of custom-designed adiabatic potentials for 
cold àtoms based on multimode radio frequency radiation in combination with static 
inhomogeneous magnètic fields. For example, the use of radio frequency combs gives 
rise to periòdic potentials acting as gratings for cold àtoms. In strong magnètic field 
gradients the lattice constant can be well below 1 /im. By changing the freqüències of 
the comb in time the gratings can easily be propagated in space, which may prové useful 
for Bragg scattering atòmic matter waves. Furthermore, almost arbitrarily shaped 
potential are possible such as disordered potentials on a scale of several 100 nm or 
lattices with a spatially varying lattice constant. The potentials can be made state 
selective and, in the case of atòmic mixtures, also species selective. This opens new 
perspectives for generating tailored quantum systems based on ultra cold single àtoms 
or degenerate atòmic and molecular quantum gases. 
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1. Introduction 

The availability of matter in a coherent state has renewed the interest in Bragg scattering 
with massive partides. Various types of gratings have been developped based on 
arrangements of two nearly degenerate làser beams pQ |2] or microfabricated structures 

These techniques have drawbacks however: Optical lattices are not flexible in the 
sense that it is difhcult to implement local variations or disorder, their geometry is frozen 
by the incident radiation of a few làser beams, and their lattice constant is limited by 
the wavelength of the employed làsers. On the other hand, microstructured gratings are 
unalterable by their design, and they are generally two-dimensional. 

In magnètic traps incident microwave or radio frequency radiation can have a strong 
impact on the trapping potential. This is due to the fact that the potential depends 
on the magnètic substate of the trapped àtoms. The energy of this substate can be 
manipulated by admixing other substates via resonant radio frequency radiation. In 
inhomogeneous magnètic fields this coupling is local and leads, within a dressed states 
picture, to avoided level crossings jl] . Àtoms moving across an area where the coupling 
is strong, follow adiabatic potentials by avoiding the crossings. 

This feature is e.g. used for evaporative cooling, where trapped magnètic states 
are coupled to untrapped states at certain distances from the trap center |ü| and 
for certain schemes of output coupling atom làser beams from trapped Bose-Einstein 
condensates [7j . Atòmic traps based on avoided crossings have been proposed by Zobay 
et al. [HI El and pioneered by Perrin et al. ^0] • Meanwhile double-well potentials have 
been built using this approach [llj . and resonators for àtoms have been demonstrated 

H2|. 

In this paper we propose to generate periòdic potentials for àtoms by combining 
appropriate magnètic field gradients with the application of multimode radio frequency 
radiation. As compared to optical or microfabricated gratings, gratings based on radio 
frequency combs possess a number of advantages. They work for àtoms with low- and 
high-field seeking magnètic moments with a phase of the grating, which depends on the 
atòmic infernal state. Furthermore, the technique allows for an almost complete control 
over the spatial shape and the temporal evolution of the potentials. In particular, 
irregular patterns may be formed, and very small structures, only limited by the size of 
the technically feasible magnètic field gradient, may give rise to gratings with very small 
lattice constants. With properly shaped magnètic fields almost arbitrary 3D potential 
geometries are imaginable. It is even possible to design time-dependent shapes, e.g. to 
shift, merge or split individual potential sites. The range of possible manipulations and 
applications seems inexhaustible. 

We begin this proposal with a quantitative description of the bàsic mechanism 
in Sec. 121 where we suggest a very simple approximation for the adiabatic potentials, 
which holds when the frequency components are not too closely spaced. To demonstrate 
the versatility of the radio frequency technique, we discuss two examples. In Sec. El 
we focus on the case of a periòdic potential (or grating). We estimate the smallest 
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periodicity, which can be realized using radio frequency combs, discuss how to generate 
moving gratings in order to perform Bragg scattering experiments and how to use the 
phenomenon of Bloch oscillations to probe such a grating. 

The second example, discussed in Sec. El is a way of modifying the shape of a 
magnètic trap using microwave radiation. In ordinary magnètic traps, atòmic spin 
states with identical Zeeman shifts experience the same confming force regardless of 
the mass of the trapped partides. The possibility to design magnètic potentials with 
arbitrary shapes provides a useful handle to selectively manipulate mixtures of different 
species. 



2. Simple theoretical description 

We consider an atòmic ground-state with hyperfme structure having the total spin F. 
The magnètic sublevels m F are split in an external magnètic field B by an amount 
^bQf^fB , where gp is the atòmic g-factor of the hyperfme level. An irradiated linearly 
polarized radio frequency B r /cos(u;t) couples the sublevels \F, m F ) <-> \F,m' F ) with 
m' F = mjf ± 1, wherever it is close to resonance, provi ded the orientations of the radio 
frequency and the magnètic fields are orthogonal. The coupling strength is given by the 
Rabi frequency [ü] 



Í2 = ^ |B r/ x è B \ y F(F + 1) - m F m' F , (1) 

where èg is the orientation of the local static magnètic field. 

Alternatively, a microwave frequency can be used to couple different hyperfine levels. 
However, the essence of the effect can be studied in the coupled system ||, §) <-> ||, — |), 
on which we will focus in most parts of this paper. A gener alization to multilevel systems 
F > ~ is straightforward. The dressed states Hamiltonian of our two-level system is a 
2x2 matrix, 

XT(y\ - ( \^B9FB{Z) - \hu \W, \ ( 

1 } ~ \ \m -^ B g F B(z) + lhu ) ■ {Z) 

For simplicity we have also assumed a one-dimensional geometry, B = B(z), but it can 
easily be generalized to three dimensions. The eigenvalues of H are 



1 



E ± (z) = ±-\M + \Mb9fB(z) - tkü] z ■ (3) 
Far enough from resonance, hfl «C \hb9fB(z) — hu\, we get 

B±W » ±\ IwfB( ,) - M ± 4bfl ^_ M , (4) 

where the second term can be interpreted as dynamic Stark shift of the energy levels. 

To illustrate the action of the radio frequency, we calculate the potential energy 
and the dressed states for 6 Li àtoms. For simplicity, we assume a 1D linear magnètic 
field gradient B(z) = zb. Fig.Ufa) shows the radio frequency coupling and Fig. [Tfb) the 
dressed states for two magnètic substates coupled by a single radio frequency component. 
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Figure 1. (a) Potential energies for a F — h level scheme with the g- factor of 
gp = — | (as is the case for the 6 Li ground state 2 Si/ 2 )- A single radio frequency 
component (arrow) couples the substates tuf = ±4- Here ò = 200 G/cm and 
lo = 2ir x 5 kHz. (b) Uncoupled dressed states (dotted line), coupled dressed states 
(sòlid), and dynamic Stark shifts (dash-dotted) calculated in the off-resonant field 
approximation. The Rabi frequency is = 2t: x 700 Hz. 



If we irradiate several frequency components u n , labelled in ascending order by 
n — 1,2, .., the problem gets more complicated, because every component has its own 
dressed states basis. However for sufficiently low Rabi freqüències and large frequency 
separation of the components the dynamics is essentially governed by the component 
which is closest to resonance. A good approximation holding for homogeneous magnètic 
fields is therefore to consider only one nearly resonant frequency and disregard the 
others. To handle the case of magnètic field gradients, we introduce a local frequency, 
which at a given location z is the frequency component closest to resonance, lü(z) = c<j n ( 2 ), 
where the component n = n(z) is chosen such that \fiBgFB(z) — %üJ n i z -)\ is smallest. 
In other words, along the field gradient we switch between different dressed state 
representations; and in order to make this switching as smooth as possible, it is done 
at those points where a blue-detuned frequency u n+ i gets closer to resonance than a 
red-detuned frequency u> n . Fig. |2fa) reproduces the energies locally dressed with the 
radio frequency field u n . 

The adiabatic potentials V a d are now obtained by considering that the couplings are 
strong enough to yield Landau-Zener transition probabilities close to unity (see below). 
Formally, this is done using the following procedure: 



EM) =F 



%ÜJ. 



n i z ) 



n{z)-l 



T E (5) 



k=l 



This is shown in Fig. |2fb). 

The adiabatic potentials in Fig. |2£b) show small discontinuities at those places, 
where the frequency component which are nearest to resonance changes. This is an 
artefact due to our neglecting of all non- resonant components. The error made can be 
reduced by considering their combined Stark shifts 

h 2 n 2 
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Figure 2. (a) Local dressed states (sòlid line). The parameters are the same as in 
Fig-H but with three freqüències irradiated at ui n = 2tt x 2,4,8 kHz. The black dots 
show the locations where the freqüències are resonant. The dash-dotted line denotes the 
frequency component uj n , which is closest to resonance. (b) Corresponding adiabatic 
potentials calculated with nearly resonant freqüències, only. (c) Dynamic Stark shift 
corrections from those radio frequency fields, which are not closest to resonance. (d) 
Adiabatic potentials with the Stark shift corrections of Fig. (c). 



in the diagonal components of the Hamiltonian 

H^z) = H„ n (z) + L n (z) ( - 1 . (7) 

The corrected eigenfrequencies are thus 

E±(z) = ± l -\jh 2 W + [fi B g F B(z) -hcü + 2L n (z)} 2 . (8) 

Fig. |2(c) visualizes the Stark shifts produced by those frequency components, which 
are not the closest to resonance. Recalculating the adiabatic potentials (J3j) with the 
corrected eigenfrequencies we obtain potentials without discontinuities as shown in 
Fig.Hd). 

A major concern is the dependence of the coupling strength on the relative 
orientation of the radio frequency polarization and the magnètic field. Any variation of 
the field vector orientation over the region, where the àtoms are trapped, introduces a 
position-dependence of the Rabi frequency. This problem can be avoided by chosing a 
sufficiently large magnètic offset field, which defines a spatially invariant quantization 
axis, and by arranging for an orthogonal polarization vector. On the other hand, this 
feature may be technically exploited as an additional handle for shaping the adiabatic 
potentials [TT] . 

The above procedure can easily be extended to an arbitrary number of frequency 
components, to higher-dimensional magnètic field shapes B(r), or to multilevel systems 
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F > |. Fig. [3^a,b) visualizes the case of a two-dimensional quadrupolar magnètic field 
|B(x, z)\ = J (zV Z B) 2 + (xV x B) 2 with three freqüències irradiated. Fig.Efc) shows the 
adiabatic potentials for five Zeeman substates, as in the case of the 7 Li ground state 
hyperfine level F = 2. Finally, Fig. Hlfd) represents the adiabatic potentials for the case 
of a frequency comb u n = 2tt x (3 + 1.5n) kHz. Note that gravity is easily taken into 
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Figure 3. (a) 2D quadrupolar magnètic field with gradients b z = 200 G/cm and 
b x — 140 G/cm. (b) Corresponding adiabatic potentials for the 6 Li two-level system 
F — \. The Rabi- and radio freqüències are íï = 2ir x 500 Hz and uj n = 2ir x 2, 4, 7 kHz, 
respectively. (c) Adiabatic potentials for a system with five Zeeman substates, F = 2. 
The <?-factor is set to gp = 5, as is the case for the 7 Li ground state 2 5*1/2. The 
Rabi- and radio freqüències are as in Fig. (b) . (d) Adiabatic potentials for a frequency 
comb and the 6 Li two-level system. In this case the Rabi frequency is chosen to be 
fi = 2tt x 400 Hz. 



account by calculating the resonant freqüències as explained above and supplementing 
the Hamiltonian ((7|) with the potential energy E pot = mgz. 



3. Radio frequency lattices 

We have seen above that by irradiating a frequency comb we obtain a periòdic potential 
[cf. Fig. EHd)]. The question arises how small the lattice constant can be made. The 
periodicity depends on the spacing of the frequency components and on the magnètic 
field gradient. The grating is characterized by the lattice constant 

_ 2h(u n+ i — Lü n ) , . 

HB9Fb 

and the potential modulation depth 

v ad = n(^±^-n) . (io) 
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To obtain small lattice constant, we chose large gradients. For the following, we set the 
premise that the maximum feasible gradient is b = 200 G/cm, which can be achieved 
with current technologies over reasonably large spatial ranges. 

A typical experiment with gratings consists in Bragg-scattering moving àtoms. 
Àtoms with a de Broglie wavelength which matches the lattice periodicity are Bragg- 
reflected. The corresponding atòmic velocity is 

md 

The adiabaticity criterion sets now an upper limit for the velocity quantified by the 
probability for Landau-Zener transitions, P^z = 1 — e^ hn2 ^ dt ^ B9FB( - z ^ with z = vt. 
The àtoms follow the adiabatic potentials if Plz — > 1- This implies híl 2 ^> fisgFvb. 
Consequently, in order to ensure Landau-Zener transitions for velocities as high as the 
Bragg velocity, the frequency comb spacing and the Rabi frequency should be chosen 
such that 



11^ 



mfl 2 d 3> HBQFb ■ 
The shaded zone in Fig. 0] depicts the adiabatic regime. 



(12) 



i 

Ci 




Figure 4. The shaded area shows the adiabatic regime for íl and d with a field 
gradient of 6 = 200 G/cm. The zone at the left of the sòlid line corresponds to 
gratings shallower than one recoil energy. The dash-dotted line delimits the interband 
tvmneling regime from the adiabatic regime (hatched area) as explained in the text. 



The interplay of the parameters íl and d also governs the potential modulation 
depth V a d through the Eqs. (jïïj) and (JTDJ). This relationship may be illustrated by 
estimating the smallest lattice constant, for which the lattice potential is still deeper 
than one recoil energy, E r = h 2 /8md 2 , as a function of íl. The sòlid line in Fig. 0] 
reproduces the curve d = d(Q) for which V a d = E r . We find that for a Rabi frequency 
of 3 kHz the lattice constant must exceed d ~ 0.7 /xm. 

The lattice constant can be reduced using larger gradients. Lattice structures on a 
length scale of a few 100 nm, which are certainly feasible, could be of interest for the 
generation of disordered potentials with the aim of studying atòmic Bose and Anderson 
glasses |13j . The most interesting regime, where the length scale of the disorder is 
smaller than the healing length of a Bose-Einstein condensate, is dimcult to reach with 
optical speckle patterns [HJ [T3]. In contrast the regime is quite accessible to radio 
frequency lattices. 
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3.1. Moving gratings 

We now want to produce an adiabatic grating which propagates in a constant magnètic 
field gradient. This necessitates a time- dependent frequency comb, ui n (t). After a 
time t n the frequency drift has covered the separation of the frequency components 
u n +i(t n ) = co> n (0), and we recover the same spectrum. The frequency separation 
determines the lattice constant through Eq. (jüj). The time t n rules the propagation 
velocity of the lattice via 

v = d/t» . (13) 
Therefore we may repeat the ramp 

u n (t) = u n (0) + {u n+1 - u n ) í" 1 (t mod t n ) . (14) 

The Fourier transform of the instantaneous power spectrum of the radio frequency, 
P{u,t) = E„í[^ - u n (t)}, is 

[/(r,í) = [/ ^ S inK(í)r]. (15) 

n 

The inverse Fourier transform of the function U(t, t), shown in Fig. Efa), is the envelope 
of the frequency comb time evolution. To retrieve the instantaneous spectrum, we build 
the invers Fourier transforms over consecutive time intervals as shown in Fig. Efb). 
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Figure 5. (a) Overall Fourier transform of the signal of 4 radio frequency components 
linearly scanned across a range corresponding to their frequency splitting. (b) Stepwise 
Fourier transform of the same signal. (c) Time-evolution of corresponding adiabatic 
potentials. The gradient is again b = 200 G/cm. Additionally, a magnètic offset field 
of Bo = 0.1 G has been applied to avoid Majorana spin flip transition at the center of 
the quadrupole trap. The Rabi frequency is Q = 2ir x 15 kHz. (d) Stepwise Fourier 
transform of a signal consisting of three frequency components, the amplitudes and 
phases of which vary differently in time. 

In the example of Fig. the chosen time base is t n = 2 ms. In a magnètic field 
gradient b = 200 G/cm, with a frequency comb chosen as uj n (0) = 2im x 100 kHz, the 
lattice constant is d — 1.7 /im. Thus the propagation velocity is v — 1.8 mm/s. 
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Multifrequency signals can be generated by a variety of techniques. The most 
convenient is probably to use arbitrary waveform generators programmed with the 
Fourier transform of the desired frequency spectrum. In order to generate moving radio 
frequency lattices, a whole period t n has to be programmed. In our example, for radio 
freqüències up to uj n Pd 2ix x 500 kHz the number of oscillations within t n is 1000. If 
we count at least 8 points to resolve an oscillation, we need to program waveforms with 
8000 points, which is feasible with state-of-the-art arbitrary waveform generators. 

3.2. Bragg reflection and Bloch oscillations 

Àtoms are Bragg-reflected at a lattice if their de Broglie wavelength is commensurable 
with the lattice constant, i.e. their velocity is given by Eq. fil 1)1 . Alternatively, for resting 
àtoms to be accelerated by a moving grating, the propagation velocity of the grating 
Eq. ([Tïïjl must coincide with the Bragg velocity. For the above parameters, we expect 
Vbrg = 1-2 mm/s. For our example, in order to obtain an interaction between the àtoms 
and the lattice, we may tune the propagation velocity by varying the time period t n or 
modifying the lattice constant d. In contrast to the propagation velocity, which increases 
with the lattice constant, the resonant Bragg velocity is inversely proportional to it. 
Hence, there can only be one field gradient satisfying the Bragg condition. Equating 
the two velocities, v = Vb rg , we derive 



which gives bb rg ~ 210 G/cm for the present conditions. 

Thus for anisotropic traps, for which the secular freqüències differ for different axis, 
a pronounced directional selectivity of the atòmic interaction with the moving lattice 
can be expected. The directions in which àtoms are accelerated are only those, where 
the Bragg condition is satisfied. 

Fig. Efc) shows a time series of the adiabatic potentials produced by the frequency 
comb specified above. When the frequency components are ramped, the potential wells 
at both sides of the trap origin migrate in opposite directions. This feature can be useful 
for realizing atom interferometers. 

The technique is obviously not limited to uniform gratings propagating with 
constant speed. As an example that arbitrary time signals may be programmed Fig.Efd) 
shows the sequential Fourier transform of a waveform consistent of three frequency 
components, which are differently modulated in amplitude and phase. 

A powerful tool to probe and characterize atom optical gratings are Bloch 
oscillations. This phenomenon occurs when a constant force, F = ma, is applied to cold 
àtoms interacting with the grating. The force can be gravity, or it can be generated by 
an accelerated motion of the grating E] • For Bloch oscillations to occur in lowest 
band, an adiabaticity criterion similar to the Landau-Zener criterion ()12|) follows from 
the Bloch model: Interband tunneling at the edge of a Brillouin zone is prevented if the 




(16) 
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rate of change of potential energy due to acceleration is smaller than the size V a d of the 
band gap, i.e. V^ d 3> ^ (hmaz)\ z=d ^ 2 . This leads immediately to 

ll^L^rnad. (17) 
4 E r 

Inserting the recoil energy and the expressions (jïïj) and (jlUp. we obtain an inequality for 
Q and d. This inequality is depicted in Fig. 0] as the hatched zone below the dash-dotted 
curve for the case that the acceleration is due to gravity a = 9.81 m/s. 

4. Potential shaping 

Another application of adiabatic potentials is the realization of non-harmonic potentials 
via the irradiation of appropriate radio freqüències and microwaves. Two-dimensional 
trapping potential based on this principle have recently been proposed [SJ |Ü]. The 
possibility to selectively influence the confinement strength of àtoms trapped in specific 
Zeeman substates or being of different species is particularly useful for studies of their 
interaction. 

An example where this plays a role is a mixture of the high-field seeking states 
6 Li ||, |) and 87 Rb |2,2) For this case, sympathetic cooling of Li via Rb to very 
low temperatures (such as those needed for reaching the critical temperature for BCS- 
pairing) requires a sensitive matching of the respective gas densities: This guarantees 
that the Li Fermi temperature is lower than the critical temperature for Bose-Einstein 
condensation of Rb and reduces fermion-hole heating as far as possible (201 • This 
respective density-matching is now possible by controlling the potentials' shape via 
microwave or radio freqüències. 

For this example, the following scenario could be used to transfer àtoms into a 
potential with a flat bottom in which an atòmic cloud adopts a nearly homogeneous 
density distribution. One starts producing an atòmic cloud, e.g. of 87 Rb àtoms in the 
|2, 2) state. With a properly designed adiabatic sweep the majority of the àtoms can be 
transferred to |2,0), and a microwave frequency resonant to the |1,— 1) state at some 
distance from the trap center is irradiated [cf. Fig. Efa)] 21 . The àtoms are then 
trapped in an adiabatic potential with a flat bottom [cf. Fig. Efb)]. Note that there is 
a second trapping potential for those |2, 2) àtoms that were initially in the center. 

5. Conclusion 

To conclude we believe that radio frequency combs in conjunction with magnètic field 
are a versatile tool for manipulating the atòmic motion on small length scales. The 
gratings are very stable and flexible. In fact the lattice constant is only limited by 
the magnètic field gradient, which is technically feasible. In our examples gradients of 
200 G/cm have been used, but gradients of several 1000 G/cm can certainly be realized 
in environments, which are compatible with cold atom òptics. Potential applications of 
time- or position-dependent radio frequency combs may include Bragg velocity filters, 
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Figure 6. (a) Hyperfine and Zeeman structure of the 87 Rb ground state. States 
|2,0) and |1, — 1) are coupled by a microwave at 6.8 GHz. The Rabi frequency is 
íi = 2ir x 600 kHz, the gradient is b = 200 G/cm, and we assume a magnètic offset field 
of Bq = 1 G. The microwave is tuned 2 MHz to the red of the field-free |1, —1) «-» |2, 0} 
resonance. (b) Dressed states give rise to two adiabatic potentials. The upper potential 
is flat around the center, the lower one has a depth of about 50 fiK for the chosen set 
of parameters. 



Bragg interferometers, quasi-random potentials with disorder occurring on a very small 
length scale, and other atom optical elements. 

In a variety of cases, such as in two-species experiments, it is important to control 
the atòmic densities without modifying the magnètic field gradients. As an example, we 
have discussed the 6 Li - 87 Rb mixture, where the feature of deforming magnètic trapping 
potentials by radio frequency radiation in a controlled manner is crucial in order for 
sympathtic cooling to be efficient down to very low temperatuires. This approach may 
thus constitute an alternative to using bichromatic optical dipole traps suggested in 
Ref. 
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